Surface Bloch waves in metamaterial and metal-dielectric superlattices 
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We study the properties of electromagnetic Bloch waves in semi-infinite periodic structures created 
by alternating metamaterial and dielectric layers. We derive and analyze the dispersion relations 
in the long-wavelength limit for both TE- and TM-polarized surface Bloch modes, for magnetic 
metamaterials with negative refraction and metal-dielectric plasmonic superlattices. We reveal that 
in the subwavelength regime the bulk modes are characterized by three different refractive indices 
("tri-refringence") while the surface modes can propagate parallel to the Bloch wavevector and along 
the interface between superlattice and semi-infinite dielectric. 



Subwavelength confinement and manipulation of light 
in nanoscale optical structures is of a great importance 
for applications in imaging, sensing, high-resolution 
lithography, all-optical signal processing, photonic fun- 
nels and superlenses However, the diffrac- 
tion limit prevents the light confinement on the scales 
smaller than a half of the wavelength. Conventional tech- 
niques for the resolution improvement are based on the 
use of high- index materials [7| , but the further progress 
in this field requires a design of artificial materials, of- 
ten termed as metamaterials, with higher refractive in- 
dices d, Plasmonic macrostructure composites that 
operate in the regime of coupled waveguide modes, have 
been recently proposed for a variety of future imaging 
systems and beam steering [Io|, US [H, HH . At the same 
time, the absence of natural magnetism at optical and in- 
frared frequencies requires again the fabrication of nanos- 
tructured metamaterials in order to achieve required dis- 
persion properties for negative refraction 0, [l5|, [H, [ljj ■ 
A typical width of a single layer in such metamateri- 
als is of the order of 10 nm, that is much smaller than 
the free-space wavelength of optical or infrared radia- 
tion. Thus, it is usually assumed that optical proper- 
ties of those multilayered composites can be described 
by an effective- medium approximation [l8| , although the 
typical field variation length may be less than the free- 
space wavelength On the other hand, many stud- 
ies of negative-index materials that use the original de- 
sign based on simultaneously negative dielectric permit- 
tivity and magnetic permeability assume that the ho- 
mogenization limit for metamaterials is fulfilled, either 
for isotropic [H US HI or 'indefinite' [H models. In 
particular, the effective-medium approximation has been 
used to study surface waves at the interfaces between 
metamaterial and surrounding media that may lead to 
the subdiffraction imaging. 

In this Letter, we study the electromagnetic waves 
propagating in two-dimensional semi-infinite periodic 
structure created by alternating layers of metamaterial 
and conventional dielectric and semi-infinite homoge- 
neous dielectric placed normally to the interfaces within 
the multilayered media [see the geometry in Fig.[lja)]. A 
unit cell of this layered structure is formed by two layers: 



one with the thickness oi, dielectric permittivity s\, and 
magnetic permeability fix, and the other one, with the 
parameters a2, £2, and U2, respectively. The layered me- 
dia has the period d = a\ + a 2 . Such a structure supports 
the propagation of surface waves in the direction paral- 
lel to the Bloch wavevector m the layered medium. 
Surface waves that propagate along the interface between 
the layered dielectric media and a semi-infinite dielectric 
cladding placed parallel to the interfaces within the mul- 
tilayered medium (i.e. with the wave vector fcj_||x per- 
pendicular to the Bloch wavevector) were studied more 
than thirty years ago in Refs. [25[ . Although many 
authors have investigated subsequently the same geom- 
etry, the surface waves that propagate parallel to the 
Bloch wavevector, have not been analyzed yet, to the 
best of our knowledge. We term those waves as surface 
Bloch waves. 

There exist surface modes supported by each individ- 
ual interface within the layered media, but the coupling 
(weak or strong) of these modes gives birth to new col- 
lective modes that propagate in the same direction, i.e. 
perpendicular to the Bloch wavevector. In the regime of 
strong coupling, when the unit cell thickness d is much 
smaller than the free-space wavelength A, the optical 
properties of a multilayered structure can be described 
by effective permittivity (l8j £ and effective permeabil- 
ity jl. Generally speaking, these quantities are diagonal 
tensors with two equal components in the directions per- 
pendicular to the Bloch wavevector: £ = diag(e j_ , e± , £ | ) ; 

fi = diag( ) tt ± ,AU,M||)- 

Assuming that the electric field of the s-polarized (TE) 
mode and the magnetic field of the p-polarized (TM) 
mode in the form: ~ Ak B (z) exp (iksz + ikj_x — iujt)y, 
the dispersion of the Bloch modes in the infinite layered 
media is described by the well-known equation for bulk 
Bloch waves that comes from the eigenvalue problem for 
two-layer transfer matrix [24], [2f| 
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where fci j2 = (£1,2/^1,2 — q'jj 1 ^ 2 - We note that all spa- 
tial dimensions are normalized to c/lv, and wavenumbers 
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FIG. 1: (a) Geometry of the problem, (b) Equifrequency 
contours for uj/uj p = 0.3 obtained from the dispersion ([2]). 



to lojc. The polarization dependent coefficients a s , p are 
defined as a s = /^Mi/^i/^j and a p = &2£i/fci£2- In the 
subwavelength regime, when |fciai| <C 1, and |fc2a2 1 <C 1, 
after the Taylor expansion in Eq. (JTJ) , we obtain the dis- 
persion relations similar to those for the light propagation 
in uniaxial crystals, for the s- and p-polarized waves 
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where n = ^e± fi± is the refractive index of the or- 
dinary wave, n s e = ^/£± M|| and r\P e = ^/sfJiT are the re- 
fractive indices of s- and p-polarized extraordinary waves, 
respectively. The effective material parameters appear as 
a result of the Taylor expansion rather than an averaging 
procedure: 



£± = (1 - S) £i + Se 2 ; £|| 
Mx = (1 +<fy*2! M|| 
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{1-8) e a + fei' 

Mi M2 
(1 - 6) fi 2 + Jmi ' 
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where 6 = (12/d, < <5 < 1 describes the degree of 
anisotropy. The structure becomes isotropic at S = 0, 
and (5 = 1, and anysotropic between this values. 

Equation @ describe the dispersion relations of bulk 
waves for infinite layered media in the long-wavelength 
approximation. For nonmagnetic media, and under the 
additional approximation \q B d\ <C 1, these e quat ions re- 
duce to the equations obtained in Refs. [241. |25||. where 
it was shown that a periodic stratified medium is for- 
mally analogous to an uniaxial crystal, yielding for s- 
and p-polarizations, respectively: q 2 B + q\ = e±, and 
q%/sx + Qj_/ £ \\ = Although these equations look very 
attractive (because they are formally analogous to the 
well known Fresnel equations), it is worth noting that 
imposing \q B d\ <C 1 leads to loosing zonal structure and 
band gaps typical for Bloch waves. Under the same ap- 
proximation, inclusion of higher order terms in the Taylor 
expansion revealed nonlocal effects (spatial dispersion) 
in nanolayered metamaterials (l9j . On the other hand, 
inclusion of magnetic metamaterials gives rise to the ap- 
pearance of the extraordinary wave even in the case of 
s-polarization, if the formal analogy with uniaxial crys- 
tals has been adopted. There are two ordinary s- and 
p-polarized modes that obey the same dispersion law: 
sin 2 qBd/2 = e±/i±(d/2) 2 and two extraordinary s- and 
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FIG. 2: Effective dielectric and magnetic parameters of the 
metamaterial superlattices (S = 0.95). 



p-polarized modes with different dispersions: q\ = £j_/X|| 
and q\ = £||M-Lj respectively. Thus, for the same fre- 
quency, we now have three different refractive indices, 
and when \qsd\ -C 1, Eqs. @ represent two ellipses in 
(<lBi gjj-plane with one common axis (see Fig.[T](b)), in- 
stead of a circle and an ellipse which we have for uniaxial 
crystals. So, magnetic metamaterials can produce either 
negative or normal (positive) refraction in all directions, 
as well as negative refraction in one direction and nor- 
mal in the other one. We emphasize that for both s- and 
p-polarizations the desired characteristics of the periodic 
stratified metamaterials can be engineered by a variation 
of the thicknesses of the constitutive layers. 

As an example, in Fig. [2] we plot the effective parame- 
ters for the case of a periodic structure with metamate- 
rial layers in vacuum (e± = fj,i = 1). Wc assume that the 
metamaterial layers have the frequency-dependent dielec- 
tric permittivity and magnetic permeability of the form: 
£ 2 = 1 - ujp/uj 2 and [i 2 = 1 — Flu 2 /(lu 2 — uj 2 ), where we 
take F = 0.56, uj 2 /iv 2 = 0.16. 

To study surface Bloch waves in long- wavelength limit, 
we assume that multilayered media is truncated, and it 
occupies a half-space x < 0, while the other half-space 
(x > 0) is occupied by homogeneous dielectric with di- 
electric constant e c . To obtain the solutions for surface 
waves, we consider in Eq. @ purely imaginary values of 
q±. Using the standard continuity conditions of the tan- 
gential field components, we than obtain the dispersion 
relation for the s-polarized surface waves, 



^ [sin 2 (q B d/2) - s^(d/2) 2 ] = -HU-J&-EC 



whereas for the p-polarized surface waves, we have 



(4) 



-!L [sin 2 (q B d/2) - e^{d/2) 2 ] = <& - e c . 

(5) 

Equations Q [5]) are obtained in the approximation 
that the thicknesses of each layer is small, however we did 
not make any assumptions about the values for the Bloch 
wavenumber. However, since we used continuous medium 
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FIG. 3: Dispersion of s- and p-polarized Bloch surface waves 
in metamaterial super lattices, for 8 — 0.95 and d — 0.1c/cj p . 
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FIG. 4: Dispersion of p-polarized Bloch surface waves in 
metal-dielectric structures, for three different values of S and 
d — O.lc/oip 



approximation for matching boundary conditions on the 
interface, it is expected that the results are valid when 
the Bloch wavelength is larger than the period of the unit 
cell (qsd < 1). More general case of arbitrary wavenum- 
ber will require numerical simulation in order to obtain 
dispersion properties of surface Bloch waves. 

Now we calculate the dispersion of the surface Bloch 
modes with the effective parameters shown in Fig. [5] for 
the waves propagating along the interface between the 
structured medium and vacuum (e c = 1). Dispersion 
of surface Bloch modes for metamaterial superlattices 
created by layers of metamaterial and conventional di- 



electric is calculated using Eqs. (|4I5[) . and it is shown 
in Fig. [3] From the dispersion equations we see that 
the p-polarized waves can exist in the frequency ranges, 
where e\\ < 0, while s-polarized waves exist for fiu < 0. 
We found that for the values of the parameter 5 close 
to 1, there are in general two branches of p-polarized 
waves and one branch of s-polarized waves. Both up- 
per and lower branches of the p-polarizcd waves are for- 
ward, while the s-polarized surface Bloch wave is back- 
ward (negative group velocity), and it exists in the range 
of frequencies 0.4 < uj/ujp < 0.603. When we decrease 
the parameter S towards the value of 1/2, the frequency 
range for the s-polarized wave existance becomes nar- 
rower, and it becomes forward (not shown), and then it 
disappears. P-polarizcd waves also disappear for values 
of 5 < 1/2, and there arc no surface waves for such values 
of 5 within the limits of our approximation. 

For nonmagnetic metamaterials or metals in the lay- 
ered structure, we use the Drude-type formula for £2(w), 
but consider /i 2 = 1. We find no s-polarized surface 
waves, but for the p-polarized surface Bloch modes we 
obtain the dispersion shown in Fig. 2J Surface Bloch 
modes exist only when the metal occupies more than 
half of the structure, i.e., 1/2 < 5 < 1. As we change 
the metal filling fraction from 5 = 1, when the frequency 
range of the mode existence is < uj/uj p < l/y/2, to 
(5=1/2, the frequency range of mode existence vanishes. 

In conclusion, we have derived and analyzed the disper- 
sion characteristics of the electromagnetic surface Bloch 
waves in metamaterial and metal-dielectric structures in 
the long- wavelength limit. We have revealed that in 
the subwavclcngth regime the bulk modes arc charac- 
terized by three different refractive indices while the sur- 
face modes can propagate parallel to the Bloch wavevec- 
tor along an interface between two superlattices. As a 
special case, we have studied the propagation along an 
interface between the superlattice and a semi-infinite di- 
electric in truncated layered structures. We note that 
the short-wavelength limit requires intensive numerical 
simulations and it will be studied elsewhere. 

The work was supported by the Australian Research 
Council. S.V. thanks Nonlinear Physics Center for hos- 
pitality and support during his stay in Canberra, and he 
acknowledges a support of the Serbian Ministry of Sci- 
ence (grant No. OI 141031). 



[1] J. B. Pendry, Phys. Rev. Lett. 85, 3966 (2000). 
[2] N. Fang, H. Lee, C. Sun, and X. Zhang, Science 308, 534 
(2005). 

[3] H. Shin and S. Fan, Phys. Rev. Lett. 96, 073907 (2006). 

[4] Z. Jacob, L. V. Alekseyev, and E. Narimanov, Opt. Ex- 
press 14, 8247 (2006). 

[5] D. R. Smith and J. B. Pendry, J. Opt. Soc. Am. B 23, 
391 (2006). 

[6] E. Ozbay, Science 311, 189 (2006). 

[7] A. W. Fang, H. Park, O. Cohen, R. Jones, M. J. Paniccia, 



and J. E. Bowers, Opt. Express 14, 9203 (2006). 
[8] A. Karalis, E. Lidorikis, M. Ibanescu, J. D. Joannopou- 

los, and M. Soljacic, Phys. Rev. Lett. 95, 063901 (2005). 
[9] A. A. Govyadinov and V. A. Podolskiy, Phys. Rev. B 73, 

155108 (2006). 

[10] A. Alu and N. Engheta, IEEE Trans. Microw. Theory 

Tech. 52, 199 (2004). 
[11] V. A. Podolskiy and E. E. Narimanov, Phys. Rev. B 71, 

201101 (2005). 

[12] R. Wangberg, J. Elser, E. E. Narimanov, and V. A. 



4 



Podolskiy, J. Opt. Soc. Am. B 23, 498 (2006). 
[13] T. J. Antosiewicz, W. M. Saj, J. Pniewski, and T. Szop- 

lik, Opt. Express 14, 3389 (2006). 
[14] I. V. Shadrivov, A. A. Sukhorukov, and Yu. S. Kivshar, 

Appl. Phys. Lett. 82, 3820 (2003). 
[15] S. Zhang, W. Fan, N. C. Panoiu, K. J. Malloy, R. M. 

Osgood, and S. R. J. Brueck, Phys. Rev. Lett. 95, 137404 

(2005). 

[16] S. Zhang, W. Fan, K. J. Malloy, S. R. J. Brueck, N. C. 
Panoiu, and R. M. Osgood, Opt. Express 13, 4922 

(2005) . 

[17] G. Dolling, C. Enkrich, M. Wegener, C. M. Soukoulis, 

and S. Linden, Science 312, 892 (2006). 
[18] A. Alu and N. Engheta, J. Opt. Soc. Am. B 23, 571 

(2006) . 



[19] J. Elser, V. A. Podolskiy, I. Salakhutdinov, and I. Avrut- 
sky, Appl. Phys. Lett. 90, 191109 (2007). 

[20] V. G. Veselago, Usp. Fiz. Nauk 92, 517 (1967) (in Rus- 
sian). 

[21] R. A. Shelby, D. R. Smith, and S. Schultz, Science 292, 
77 (2001). 

[22] I. V. Lindell, S. A. Tretyakov, K. I. Nikoskinen, and S. 

Ilvonen, Microwave and Opt. Techn. Lett. 31, 129 (2001). 
[23] D. R. Smith and D. Schurig, Phys. Rev. Lett. 90, 077405 

(2003). 

[24] P. Yeh, A. Yariv, and C. S. Hong, J. Opt. Soc. Am. 67, 
423 (1977). 

[25] A. Yariv and P. Yeh, J. Opt. Soc. Am. 67, 438 (1977). 



